Influence of the passive region on Zero Field Steps for window Josephson junctions 
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We present a numerical and analytic study of the influence of the passive region on fluxon dy- 
namics in a window junction. We examine the effect of the extension of the passive region and its 
electromagnetic characteristics, its surface inductance Lj and capacitance Ci. When the velocity in 
the passive region vi is equal to the Swihart velocity (1) a one dimensional model describes well the 
operation of the device. When vj is different from 1, the fluxon adapts its velocity to vi. In both 
cases we give simple formulas for the position of the limiting voltage of the zero field steps. Large 
values of Lj and Cj lead to different types of solutions which are analyzed. 

PACS numbers: 85.25.Cp, 84.40,74.50.+r,03.40.Kf 



I. INTRODUCTION 

The present design of low temperature superconduct- 
ing devices based on semi-conductor technology enables 
to integrate Josephson junctions into superconducting 
strip-lines to realize complex devices. A simple example 
is the window design shown in Figure 1 where a single 
rectangular junction is surrounded by a uniform passive 
region. In other systems like integrated receivers [^] , the 
local flux flow oscillator is coupled directly to a strip-line 
containing a few junctions which will realize the mixing 
of the local signal with the external signal coming from 
an antenna. 

The electrodynamic behavior of such window junctions 
is given by the phase difference between the top and bot- 
tom superconducting layers. It obeys an inhomogeneous 
sine-Gordon equation which is equivalent to a homoge- 
neous sine-Gordon equation for the junction domain cou- 
pled to a wave equation for the passive region. 
In the static case, a small passive region causes a rescal- 
ing of the Josephson length Aj into AeS > Aj |^ and 
this renormalization was seen in the experiments by the 
group of Ustinov |^ . On the other hand a large passive 
region causes the destruction of the kink due to the finite 
length of the junction [Q. 

In the dynamical case two limiting systems can be con- 
sidered, first the passive region can be present only along 
the longitudinal direction of the junction. Then one can 
assume a transverse profile that propagates in the long 
direction. Using these ideas, Lee et al derived the disper- 
sion relation for linear superconducting strip-lines ^. 
The motion of kinks in such a window junction was stud- 
ied by one of the authors using periodic boundary con- 
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ditions in x ||^. It was shown that the speed of the kink 
depends on the extension of the passive region and that 
Cerenkov resonances occur between the cavity modes and 
the kink. These show up as steps in the zero field steps 
(ZFS) in the IV characteristics. The experiments con- 
ducted in this geometry, for a fixed set of electric param- 
eters confirm these findings. 

In the second limiting case, which we will call the ID 
model, the passive region exists only at each end of the 
junction as shown in the bottom left panel of Figure 1. 
Then the passive region acts in a different way, and fixes 
the speed of the kink ||] . In |^ we also studied the 2D 
problem i.e. a homogeneous passive region surrounding 
the junction, using an adapted numerical method, the 
finite volume approximation together with soliton per- 
turbation theory. The ID model showed that interfaces 
can act as fluxon traps. Numerical results in the two di- 
mensional case indicated that the sheet inductance miss- 
match Lj I Lj between the junction and the passive region 
is the main cause of instability of the kink motion. 

Here we confirm these results, give the details of the 
characteristic IV curves and ZFS for the ID and 2D mod- 
els and compare these two situations for different geomet- 
rical and electrical parameters. We also give a complete 
derivation of the continuum model used in j|] from the 
Resistive Shunted Junction (RSJ) approximation and ex- 
plain discreteness effects observed in the solution of the 
ID model. For the 2D model we justify the fact that the 
kink adapts its velocity to the one in the passive region. 
Finally we provide insight into the kink stability for large 
sheet inductance in the passive region L/ or capacitance 
per unit surface C/. In the 2D case we give specifically 
the different IV curves which can be seen for the different 
values of the sheet inductance in the passive region Lj. 

The paper is organized as follows. After deriving the 
continuum equations for the RSJ model in section 1, we 
introduce the ID model and study its zero field steps in 
section 2. In Section 3 we study the ZFS for the 2D 
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model and compare them to the ones for the ID model. 
We give our concluding remarks in section 4. 
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FIG. 1: A view of a window Josephson junction. The bottom 
left panel shows a schematic top view. For the system shown 
on the right the linear region exists only on the left and right 
sides of the junction. 




FIG. 2: Portion of equivalent circuit model corresponding to 
the window junction displayed in Fig. 0. Crosses represent 
tunnel elements. The critical curent If is nonzero only inside 
de window 



II. 



THE MODEL 



A. RSJ model for window Josephson junction 

A simple mathematical model of the window Josephson 
junction is to describe each superconductor by an array 
of inductances L (see Fig. ||). The coupling elements 
between two adjacent nodes in each array are, a capacitor 
C, a resistance R and a Josephson current Ic |l^ . The 
Kirchoff laws at each couple of nodes (if,, h) in the bottom 
and top superconducting layers can be combined to give 



the relation expressing the conservation of currents at 
node i in the Josephson junction 



3 ■' 

and in the passive region 



0, (1) 



(2) 



where $ = — $b (resp. = vfj — vj/^)^ is the phase dif- 
ference between the two superconductors in the junction 
(resp. passive) part, the summation is applied to the 
nearest neighbors. Also note that equations (|l|) and (|^) 
are discretisations of Maxwell's equations (wave equation 
part) and Josephson constitutive equations (sinus term), 
assuming an electric field normal to the plates, a mag- 
netic field in the junction plane and perfect symmetry 
between the top and bottom superconducting layers. We 
can now obtain the model in the continuum limit, more 
suitable for analysis. 



B. Continuum limit 

The continuum version of the system (|l|)-(^) can be 
derived by introducing the following quantities per unit 
area {a?) of elementary cells of length a. 
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I, J 



I, J 



Jc 



a" a" 
We normalize the phases by the flux quantum $o. 
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and introduce the Josephson characteristic length Aj = 
a\,j 
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Notice that in this 2D problem the inductance associated 
to each cell is equal to the branch inductance Lj (resp. 
Lj) in the junction (resp. passive region). This is not the 
case for the ID model for which these inductances are in 
series, giving a total inductance proportional to the mesh 
size dx. 

We substitute relations (|) - (|) in (0) and (1) and 
obtain in the junction 



a Jc sm ( 
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and in the passive part 



= 



(6) 



(7) 



3 



After simplification of Eqs. (||)-(|7|) and using Eq. (||) we 
obtain 



CjL 



A} 
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= . 
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Now we take the limit a ^ 0, so that 



—Acj) and ■ 
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and obtain the following system of partial differential 
equations 



A0+J^sin0+^f .0 



(10) 



CiL 



= 



(11) 



We introduce the plasma frequency in the junction ojj = 
Ajwj^ , where vj'^ = CjLj is the Swihart velocity [^ . 
Similarly we define the electromagnetic wave frequency 
in the passive part uij'^ = AjW~^ , where v~'^ = CiLj. In 
the following we fix the couple (Lj, Cj) and loj and vary 
the inductance Lj and capacity C/ in the passive region. 
We normalize equations ( p^ - (|ll|) in space by Aj and 
in time by uJj^ and obtain in the junction {flj) 



sm ( 



and in the passive region 



(12) 



(13) 



where the dimensionless parameter a = 4?" a is 

V Cj 

a damping coefficient which depends on the spatial 
discretization via Aj. 



The boundary conditions of equations (^_2|) and ( |13| ) 
are of inhomogeneous Neumann type which physically 
indicates a lateral injection of current or/and an external 
magnetic field 



— \7ip ■ n = Jext on dfl 
Li 



(14) 



where n is the exterior normal. To this we add the in- 
terface conditions for the phase and its normal gradient, 
the surface current on the junction boundary dflj 



I di) d(j) 
il) = (j) and — — = — 
Li an on 



(15) 



Note that the jump condition (2nd relation in (|T^)) can 
be obtained by integrating (12) - ( |l3|) on a small surface 
overlapping the junction domain 17 j. 
In the rest of the paper we assume a rectangular window 
of length / = 10 and width w — \ embedded in a rectan- 
gular passive region of extension w' as shown in Figure 
1 . We will not consider the influence of an external mag- 
netic field and will assume the external current feed to 
be of overlap type so that the boundary conditions ( p^ ) 
become 
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where C — 2w' 4- 1 and W = 2zi;' -I- w. 



We assumed throughout the study a small damping 
a = 0.01 which is typical of under damped Josephson 
junctions. 



III. A ID WINDOW JOSEPHSON JUNCTION 

A. The model 

We first consider the simplified situation where the pas- 
sive region is present only at the two lateral ends of the 
device as shown in the bottom left panel of Fig. 1. For 
this system it is possible to derive the effective ID equa- 
tion 

C{x)tftt - dx 1 TT^'/'a; ) + e(a;)( sin^ + a(pt) ~ j{x) = 

for xe]0,C[, (17) 

where ip is the phase averaged in y across the junction 
and 7(2;) = L(i.)s where s is the junction area. 
In equation ( p7| ) C, L and e are discontinuous functions 
of X defined by 



iC{x),L{x),e{x)) = 



(1,1,1) 



if X (e]w',w' + l[ 



(Ci.Li, 0) if X e]0, w'[U]w' + l,C[, 



The boundary conditions (|l6|) reduce then to 

(p,(o) = (^,(/:) = o. 



(18) 



(19) 



To compare different values of w' we normalize the 
current / by the maximum current that the junction can 
carry /max = s. To simplify the calculations we have 
assumed a constant current density 7. 

We choose the values of the inductance L/ and capacity 
Ci according to two strategies. First we assume Li Ci = 
vj"^ is constant which gives a uniform velocity in the 
passive part. In a second set of calculations we vary vi 
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by moving orthogonally to the set of hyperbolas in the 
(L/, C/) plane. 

The numerical method to obtain the (IV) curve is to 
start with a static kink like solution obtained by solv- 
ing the static problem derived from ( [l7| ) and to increase 
progressively the current to obtain a moving fluxon. Af- 
ter obtaining a stable solution in time, we compute the 
voltage in the center of the junction by 



V^—^ f\t[w' + l/2,t)At = 

^{w' + l/2,T2) - ipjw' + l/2,Ti) 
T2-T1 



(20) 



where the times Ti et T2 are given such that the voltage 
is stable. 



B. Discreteness effects 

In this section we consider a window junction with ho- 
mogeneous electric properties [Li = C/ = 1) and fix 
the width of the passive part w' — 2. We have chosen 
the number of discretization points in the x direction as 
N = 253, 503, 703 and 1003. In this case we compute the 
IV characteristic as shown in the top panel of Fig. ^. 
For small values of N, we observe resonances in the 

IV curve which decrease and disappear as the number of 
mesh points is increased from N ~ 253 to 703. These 
disappear altogether for N = 1003. To see if these reso- 
nances are due to the presence of the passive region, we 
have calculated a IV curve for a homogeneous junction 
{w' = 0) and indeed they are absent. To understand 
the mechanism of this fine structure in the IV curve, we 
have plotted on the bottom panel of FigjH the instanta- 
neous voltage (f>t as a function of x for a fixed time and 
N ~ 703, 1003 and 1403. From the comparison of the 
two panels of Fig. |^, it is clear that the radiation is re- 
sponsible for the fine structure. In fact the wave length 
of the excited radiation corresponds to a standing mode 

V = 2 mr/4: with n « 10. 

Let us now explain how such radiation is excited by the 
kink as it crosses an interface from a linear to a nonlinear 
medium. Fig. ^ shows the phase velocity = u)/k 
for the continuum sine-Gordon system (where su\<f> k, 
(f>) and the discrete system used numerically for which 
the second derivative is approximated by a three points 
difference (pxx = ('/'n+i + 0n-i — 20„). For the former 
Vtfi = a/1 -I- l/fc2 > 1 so that nonlinear modes cannot 
be excited via the long wave short wave Benney resonance 
mechanism |jl^. This is not true for the discrete system 
for which 



W0 = 



1 



/l2 



sm'^{kh/2) 



can be smaller than 1 for h and small wave lengths. 
As expected the threshold k decreases as is increased. 



Notice that the scale in k corresponds to 9 x 7r/2 < k < 
12 X tt/2. 

We then explain quantitatively the resonances in the 

IV curve shown in Fig. || by the fact that as the 
kink velocity increases, it can lock with a given cav- 
ity mode. For example N = 703 gives voltage (resp. 
velocity) steps V = 0.4457,0.44645,0.447,0.4475 (resp. 

V = 0.9935, 0.9952, 0.996, 0.997) which correspond to the 
cavity mode index n = 14, 13, 12 and 11. For N = 503 we 
observe steps at voltages V = 0.4451, 0.446, 0.4467 (resp. 
V = 0.9922,0.99426,0.9958) for which n = 11,10 and 9 
respectively. The lower values of n explain the larger 
amplitudes observed in the radiation. 

An intrinsic fine structure in the IV characteristics has 
been seen in experiments with window junctions with a 
lateral passive region. For example Thyssen et al. in |l^ 
show a slight shift in the position (and therefore limiting 
fiuxon velocity) of the IV curves for different values of w' . 
Ustinov et al Q also display IV curves which depend on 
w' and large substructures. Recent numerical work on 
a Josephson window junction with a lateral passive region 
and periodic boundary conditions along the propagation 
direction confirm these resonances due to Cerenkov radi- 
ation between a soliton travelling faster than v = 1 and 
the radiation of phase speed — + > 1. Here 
we show that when the passive region exists only in the 
propagation direction, resonances disappear at large res- 
olutions. This can be expected from the calculation of 
the emitt ed p ower of radiation for a kink as it crosses an 
interface |l6| . This quantity drops fast to zero as the 
speed approaches 1. 

In the rest of the study we have chosen N = 1403 
so that numerical resonances are eliminated from the IV 
characteristic. 



C. Zero field steps (ZFS) 

We now proceed to give an estimate of the time average 
of the voltage for the window junction. For this, we com- 
pute the velocity in the junction using the McLaughlin- 
Scott soliton perturbation equations p7[ and get for the 
kink velocity 



1 



When the current 7 increases and becomes larger than 
a, the velocity v gets close to 1. 



In the passive part the velocity vi is given by (17) 

1 



Vi 



so that the average of the voltage in time in the limit of 
large current can be approximated by 



ID 



;< Vt >-- 



Aiyj _ 27r 

~ 2w'/vi + l/v 



(21) 
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FIG. 4: Phase velocity of linear waves for the discrete model 
vs. k for three values of the step discretization h = 
0.01, 0.02 and h = 0.03 



large value of w' (typically w' — 12), the kink becomes 
unstable. This seems to be due to the fixed points that 
exist on each junction/passive region interface The 
kink can be trapped easier by one of them as w' increases 
because the driving force due to the current gets weaker. 



FIG. 3: ID dynamics. The top panel shows the current volt- 
age characteristic for w' = 2 and various number of discretiza- 
tion points. The bottom panel shows the corresponding in- 
stantaneous voltage 'fit{x) for a given time for three different 
values of the discretization A'^ — 703, 1003 and 1403. The 
current density is equal to 0.08. 



where the denominator is just the time take by the kink 
(phase jump of 2tt) to travel across the device. We now 
confirm this estimation first when the electric properties 
of the device are homogeneous and then when they are 
different in the junction and the passive region so that 
VI ^ 1. 



D. Influence of geometry 

First we consider a situation where the electric proper- 
ties are homogeneous in the whole device so that (Lj = 
Cj — 1) and change the extension of the passive region 
w'. 

In Fig. ^ we plot the (IV) characteristic curves for four 
values of w' — 0,1,2 and 3. This figure, show that the 
position of the ZFS moves towards the left when w' in- 
creases. This agrees with formula (^l|). Table (|) reports 
the limiting values of the zero field steps and the esti- 
mates from (pl|), which are in excellent agreement. For 
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FIG. 5: The ID current voltage characteristic for Li = Ci 
and various extensions of the passive region w' . 



E. Influence of the electrical parameters 

The second case tested is when the velocity in the pas- 
sive region vj 1. To simplify things we fix w' to 2 in all 
the presented runs. We show in Fig. ^ the IV character- 
istics for four values of the velocity in the passive region 
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TABLE L positions of the ZFS for the ID window model 



homogeneous (L/ 
(lit)- 



Ci = 1) together with the estimates 



TABLE II: positions of the ZFS for the ID window model 
together with their corresponding estimates for w' — 2. 





LiCi 


VI 


VlD 


VnurnlD 


w' 


VlD 


VnumlD 




0.25 


2 


0.523 


0.521 





0.628 


0.625 




1 


1 


0.449 


0.448 


1 


0.523 


0.521 




4 


0.5 


0.349 


0.348 


2 


0.448 


0.447 




9 


0.33 


0.285 


0.284 


3 


0.392 


0.391 





F. Solutions for large Lj or Ci 



vi = 0.33, 0.5, 1 and 2. Notice that the positions of the 
ZFS of the ID effective model follow as shown in 
table Q). Fig. ^ also shows that the position of the ZFS 
moves from the right to the left as vj decreases from 2 to 
0.33. When vj goes to zero so that Lj or Ci goes to infin- 
ity, we expect the ZFS to disappear. We observe regions 
of instability for L/ >> 1 or C/ >> 1 and the zero field 
steps exist only in small intervals in the {Li,Ci) plane. 
In Fig. ^ we present a detailed numerical exploration 
of the {Li,Ci) plane. The hyperbolas corresponding to 
LjCj = w/"^ constant are shown. We indicate by the 
symbols (*), (x) and (+) the points where the ZFS ex- 
ist. We note that for these points the solution of the 
window ID problem is a kink. The region of stability is 
concentrated along the diagonal. Note that the results of 
Fig. |8| are based only on the numerics. 
When C/ or Lj become large we obtain particular types 
of solutions which we now discuss. 
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Instabilities of the kink motion occur when the exten- 
sion of the passive region w' is large, typically w' ~ 12. 
In that case the current density becomes small and can- 
not drive the kink out of the potential well created by the 
interface Q. Other instability factors are the electrical 
parameters (L/,C/). We have shown in |9[ |l^ that the 
electrical parameters act differently from the geometry. 
For Ci >> 1 the solution is static while it is dynamic for 
L/ >> 1. In this latter case, kink motion is possible only 
in the junction. 

We now consider the disappearance of the Zero field steps 
in the IV curves. This occurs for intermediate values of 
C/ and Li and gives rise to a large voltage. 
Let us recall the equations of the problem 



0.6 



g2. 

9t2 



52 



a— w < X < w + I , 
ot 



C 



— J < X < w' and 



w 



(22) 

-I <x <C 
(23) 



together with the interface condition at x — w' ,w' + 1 

1 dip d(j> 



Li dx dx 



and homogeneous boundary conditions ^\x=o — 

dx \^—^ ^' 

In the limit of large voltages V 

-^V and -^V 

so that one can average the above equations which then 
reduce to 



dx 



^ ^ _ Q,y w' <x<w' + 1 , (24) 



FIG. 6: The ID Current voltage characteristic for w' — 2 and 
various velocity in the passive region. The velocities vi are 
indicated in the top right corner of the picture. 



1 

Li i9a;2 



7 < a; < w' and w' + I < x < C , 



(25) 



where the voltage V = '~fC/{al) can be derived from the 
work equation (see Appendix). 
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The time dependent part of (f>, V' is now just Vt and the 
X dependent part verifies a boundary value problem with 
'4'x\x=o,c = and the above given interface conditions. 
The solution is symmetric with respect to x = C/2 and 
we obtain for the left and right passive regions and the 
j miction 



Vt, 



V;, = -7l /^ + Vt 



(26) 



(27) 



7// ?// '2,11)^ 

^^-^-x{C~x)+^—{2-Li + —) + Vt (28) 

The above given expressions give a very good approx- 
imation of the phase computed numerically as can be 
seen in Fig. ^ which shows 'ip{x) and (l){x) for 11 suc- 
cessive values of time. The electrical parameters are 
C/ = 1 , L/ = 10 and the top panel corresponds to a 
small current density 7 = 0.1 while the bottom panel 
is for 7 = 0.25. The agreement for the latter is excel- 
lent, showing complete overlap between the approxima- 
tion (^,^,^8|) and the numerical solution. On the top 
panel the approximations give the overall value but there 
are some oscillations due to a resonance between the main 
frequency V and a subharmonic. 

It is possible to estimate the validity of the approx- 
imate solution ( [26|j27|j2^ ) by considering perturbations 
around it. This leads to the linearized equations for the 
perturbations 



9i2 



dx"^ 



- ^ cos ( 



w' < X < w' 



+ 1. 



(29) 
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Li dx^ 



< a; < and 



w 



I <x<C 
(30) 



with the interface condition at x ~ w' + I 



J_dC _ 
Li dx dx 



and boundary conditions ^|a;=o = fjU=c — 0. 

We then assume a uniform time dependence e~"*/2+*'^* 
and obtain a spectral problem which can be solved. We 
can neglect the damping because it is small. 

In the particular case presented in the top panel of Fig. 
the voltage is 1^ = 14 and the subharmonic observed 
\sV/A — 3.5. This corresponds to the mode rij = 11 in 
the junction and n/ = 7 in each lateral passive region of 
extension 2, taking into account the inductance Li = 10. 

Fig. H shows the regions of existence of ZFS in the 
{Li,Ci) plane. Notice the symmetry. For C/ >> 1 the 
solution is static while it is dynamic for Li >> 1, this 
latter parameter plays a predominant role due to the in- 
terface conditions. 




10 12 14 




10 12 14 



FIG. 7: Phase (fi versus x for eleven successive times separated 
by At = 0.11. The parameters are Ci — 1, Lj = 10, 7 — 0.1 
V — 14. (top panel) and 7 — 0.25, V = 35. (bottom panel). 
The numerical results are given in solid lines while the analytic 
expression ( p^ , p7|j2^ ) are given in dashed lines. 



IV. 2D WINDOW JUNCTION 

In this section we discuss the dynamics of the kink in 
the window junction and how to reduce it to a ID ef- 
fective problem. It is known that for a pure Josephson 
junction this reduction is possible. For example in the 
static case the behavior of a junction of width w < 2tt 
is very well approximated by the equation for the zero 
order mode |1^. For the dynamical case Eilbeck et al. 
pof found that when {I/8C <C 1) the phase is uniform 
in the y direction so that the reduction is possible. For 
a window junction the situation is different because we 
have a passive region around the junction which has elec- 
trical properties, a capacity and an inductance which fix 
the velocity vi . The comparison between 2D and ID was 
done for the two cases of homogeneous electrical prop- 
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The velocity v of the kink on the boundary of the junction 
is defined by 



on dflj . 



(34) 



and the same for the velocity vj of the linear wave on the 
boundary of the junction 



vj — on Oil J 

Vx 



(35) 



Thus from ( |3l| ) and (pl^ ) we obtain v = vi so that the 
estimate for the average voltage is then 



/111' 



27r 



2w' /vi + l/vi 



(36) 



FIG. 8: Parameter plane (Lj, Ci) showing the regions of exis- 
tence of zero field steps (ZFS) corresponding to the shuttling 
motion of a fluxon in the ID case. The velocities vi are indi- 
cated in the top right corner of the picture. 



erties in different geometries and a fixed geometry with 
different electrical properties. Before giving these results 
we discuss how the kink velocity is fixed by the presence 
of the passive region. 



A. Justification of the effective ID model 

When the velocity in the passive region w/ 7^ 1 the for- 
mula ( I21] ) is completely off for the 2D calculations and 
the reduction of the 2D window junction problem to a 
ID effective problem is impossible. 

Here the lateral passive region where waves can only 
propagate at velocity vi is controlling the kink mo- 
tion. The velocity -a free parameter for the sine- 
Gordon equation- can then adjust itself to vi so that the 
"dressed" kink can propagate in the device (junction and 
lateral passive region). The fluxon adapts its velocity to 
vi so that V = vj. 

To justify this, we take the time derivative of the first 
interface condition ( |l5| ) and obtain 



(pt = -0t on dflj , 
the second interface condition (|l^) reads 
V0 • n — Vip ■ n on dflj 



(31) 
(32) 



The numerical simulations show that it is legitimate to 
assume that the solution is a ID kink travelling in the 
passive region with the velocity vi and in the junction 
with velocity v. Then ( ^ ) becomes at the longitudinal 
boundaries 




FIG. 9: Contour of the phase for a homogeneous window 
junction with w' = 2 and / — 0.22 



B. Influence of the geometry on the zero field 
steps 

In this case the electrical properties are homogeneous 
in the whole device so that {Lj — Cj = 1). The numeri- 
cal simulations show that for different w' the solution of 
the 2D problem are y— uniform. For w' — 2 for example 
the Fig. H represent the contour of the phase in the whole 
domain at the time T = 5.024 10''. We note that the Fig. 

was normalized by max|0|. The levels lines are parallel 
in the y direction so that the solution is very close to a 
ID kink. This ID kink propagates with a velocity vj. 

The figure ( p^ ) shows the IV characteristics for the 
both ID and 2D problems. One can see the good agree- 
ment between the ID and 2D calculations. Table (HI) 



I 



dxij) = d^ip on di}j 



(33) 



reports the limiting values of the zero field step com- 
puted numerically and the corresponding estimates from 
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FIG. 10: Current voltage characteristic for a homogeneous 
window junction and various extensions of the passive region 
w' . Both ID and 2D calculations are presented for compari- 
son. 



( pl| ) in excellent agreement. 

As the numerical results show, in this case the dynamics 
of the solution in the 2D window problem can be de- 
scribed by a unidimensional fluxon and the reduction to 
a ID effective problem is possible. 



TABLE III: positions of the ZFS for a both ID and 2D ho- 
mogeneous junction together with the estimates (|2l|). 



w> 


VlD 




Vnum2D 





0.628 


0.625 


0.606 


1 


0.523 


0.522 


0.510 


2 


0.449 


0.448 


0.442 


3 


0.392 


0.391 


0.386 



C. Influence of the electrical parameters 

We have obtained the same average voltage for differ- 
ent values of Lj Ci as long as the product Lj Ci = vj'^ is 
kept constant. We observe that the position of the zero 
field step moves from the right to the left as the velocity 
in the passive region vj decreases from 2 to 0.5. We ex- 
pect the ZFS to disappear when vj goes to ie Lj or Cj 
go to infinity. 



D. 



Instabilities 



As in ID case there are two types of instabilities of the 
kink. The first is due to the geometry of the window, 
because large widths w' of the passive region (typically 
w' = 7) lead to a "stretched" kink in the junction which 



strongly radiates. It is the impossible to accelerate such 
an object and have a stable zero field step 1^. This solu- 
tion which is almost uniform in the junction will oscillate 
but not rotate and the average voltage will be zero [||. 
This situation was observed in the experimental work of 
Thyssen et al. who never found any zero field step 
for w' /w < 3. 

The second type of instability is due to the electrical 
quantities Lj and C/. For w' — 2, we have investigated 
the {Li,Ci) plane and found numerically the regions of 
instability of the kink as shown in Fig. |l^. We also 
plot three hyperbolas corresponding to vj = 0.5(*), l(x) 
and 2(-|-). The markers indicate the positions where we 
have found the ZFS. All the points for a given marker 
correspond to the same voltage. In Fig. |l^ we have 
isolated the stability region as the interior of the domain 
bounded by the solid line. We observe that the instability 
is due essentially to the inductance Lj so that the kink 
motion is stable for 0.1 < Lj < 2, independently of C/. 
For example C/ = lO"' and Lj = 1 gives rise to a ZFS. 
We can estimate the average voltage using formula (|3^), 
(V2D = 4.5 X 10~^ which is very close to the value found 



numerically Kum2£' = 4.577 x lO^"*). Thus even for 

such a large value of C/, the formula (|3^) gives a good 

agreement with the numerical calculation. 

Small values of C/ lead to a static solution as can be seen 

from equation ( [l3| ) which reduces to Aip = so that ip 

is constant. Because of the interface conditions (p is also 

constant. 

The influence of the inductance Lj has been analyzed 
in detail by calculating IV characteristics for w' — 3 and 
0.1 < I// < 3 which are plotted in Fig. |l2|. There are 
three main regions identified by the letters A B and C 
which correspond to three different dynamical behaviors. 
In region A 2 < L/ < 3, we obtain a static solution and a 
zero voltage. At this time we do not understand why the 
kink motion becomes unstable for Lj > 2 independently 
of C/. This effect is present only in two dimensions and 
therefore linked to the presence of a lateral passive re- 
gion. It may be due to the large jump in the gradient 
(j)y — {1/Lj)'tl>y along the lateral interface. 
Region B corresponds to a ZFS thus a stable kink mo- 
tion in the device. The limiting voltage is given to a very 
good approximation by formula ( |3^ ) as shown by the in- 
set of Fig. [iJ. This expression also gives a good estimate 
even for vi > 1. Region C Lj = 0.1,0.2 corresponds to 
much larger values of the voltage and a straight behavior 
typical of linear resonances in the junction cavity. Small 
values of Lj give rise to a spatially uniform phase in the 
passive region, which reacts to the phase in the junc- 
tion. The voltage for Lj ~ 0.1 corresponds to a higher 

67r/10. 



13 shows 



order cavity mode V = 2.56 ~ 67r/10. Fig. 
the corresponding phase at a given instant (top panel) 
together with the temporal behavior (bottom panel) . As 
expected the phase is almost uniform, the bounds in Fig. 
[1^ are 191.83 and 193.14. The phase increases linearly 
with time. 

Then we conclude that the region of stability of the 
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kink in the window problem is 

Atabiiity = {0.1 < i/ < 2 and C/>0.5}. (37) 
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FIG. 11: Parameter plane {Li,Ci) showing the regions of ex- 
istence of zero field steps (ZFS) corresponding to the shuttling 
motion of a fluxon for 2D case. The velocities vi are indicated 
in the top right corner of the picture. The parameter w' = 2 
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FIG. 12: IV characteristics for 0.1 < L/ < 3, C/ = 1 and 
w' = 3. The insert shows the fluxon velocity in region B as a 
function of the linear wave speed vi in the passive region. 



V. CONCLUSION 

In this paper we have investigate the electrodynam- 
ics of a uniform window junction, using the limiting case 
termed ID model where the lateral extension of the pas- 
sive region is neglected. We use the solution of the ID 
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FIG. 13: Three dimensional view of the phase at a given time 
(t=30800.) for Li = 0.1, region C in the previous picture (top 
panel). The bottom panel gives the temporal behavior of the 
phase in the middle of the junction. 



model as a guide to discuss the behavior of the full win- 
dow junction and this enables us to give an estimate of 
the average voltage for different extensions w' of the pas- 
sive region and different electrical parameters which fix 
the velocity in the passive region vi. This estimate is 
in excellent agreement with the numerical results even 
when vj ^ I. In this situation we show that the fluxon 
adapts its velocity to vi. For both cases v/ = 1 and 
Vj ^ 1 the numerical results show that the fluxon prop- 
agates perpendicularly to the direction of the junction. 
We have doubled the jmiction length to / = 20 and find 
the same estimates for the limiting voltage. In the 2D 
window junction, the translational symmetry is broken 
so that we obtain a different result from the case where 
only a lateral passive region exists. 

We have also studied numerically the stability of the 
zero field step in the window junction and show that when 
the width of the passive part becomes large the fluxon 
becomes distorted and gives to another type of solution 
which is radial. In this case we have not found any ZFS. 

This study shows that the kink can travel into the pas- 
sive region even when the impedance is not adapted as 
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long as the mismatch is in the capacitance and not the 
inductance. In other words the hmit e — *■ is not a sin- 
gular limit. On the other hand if there is an inductance 
mismatch then this will most likely cause the break-up 
of kink shuttling leading to the disappearance of the zero 
field step. The stability of the fluxon depends essentially 
on the inductance Li. This result, we believe, can lead to 
improving Josephson devices and their coupling to micro- 
strip lines. 
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Integrating by parts the second term in the left-hand side 
and using the boundary conditions we obtain 



d_ 

di 



C{X) 2 



1 



2L{x) 
c 



'y{x)ipt dx — a ^{x)'^t da; ■ 
Jo 



Recalling the sine- Gordon- wave Hamiltonian 



H = 



C{x) ^ 



1 



2L{x) 



(p^ +e(a:) (1 - cos(^) 



da; 
(Al) 



The computation of the linear part of the IV curve is 
based on the conservation of the ene rgy of the system 
that is given by sG-wave Hamiltonian (Al) from which it 
can be shown that H satisfies the power balance equation 



APPENDIX A: DERIVATION OF THE LINEAR 
PART OF THE IV CURVE 



AH , '•"''+' 

— — = 7 / Ptax- a 

In the stationary regime the average dH/dt = , so 



ifl Ax 



(A2) 



To determine the Hamiltonian associated with the ID- 
window junction equation we multiply Eq. ( |l7| ) by dtp/ dt 
and integrate between and C to obtain 



d_ 
di 

r 



L{x) 



Lp^ + e{x) (1 — cos Lp) Ax 



(fit Ax 



j{x)ipt Ax - 



a I e{x)(pt dx , 
/o 



J I PtAx 



rw'+l \ 

/ ^?da;)=0. 



Since the voltage V is the temporal avera|;e of ipt, we 
obtain in the high voltage limit (^ipl) = {ipt) , so that the 
IV curve is given by 
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